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WAS THE UNIVERSE REIONIZED BY MASSIVE METAL FREE STARS? 
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ABSTRACT 

The WMAP satellite has measured a large optical depth to electron scattering after cosmological 
recombination of t^s = 0.17 ± 0.04, implying significant reionization of the primordial gas only 200 
million years after the big bang. However, the most recent overlap of intergalactic HII regions must 
have occurred at z < 9 based on the Lya forest constraint on the thermal history of the intergalactic 
medium. Here we argue that a first generation of metal-free stars with a heavy (rather than Salpeter) 
mass function is therefore required to account for much of the inferred optical depth. This conclusion 
holds if feedback regulates star formation in early dwarf galaxies as observed in present-day dwarfs. 

Subject headings: cosmology: theory - early universe - intergalactic medium - Stars: formation 



1. INTRODUCTION 

The first stars in the universe (Pop-Ill) formed out of 
metal-free gas, relic from the big bang. The bottom- 
up hierarchy of structure formation in Cold Dark Mat- 
ter (CDM) models implies that stars started to form in 
CDM halos just above the cosmological Jeans mass of 
~ 1O^M0 at redshifts as high as z ~ 20 — 30 (see review by 
Barkana & Loeb 2001 and references therein). Simulations 
of metal-free star formation indicate that the first stars had 
high masses (M > lOOM©) since gas cooling by molecu- 
lar hydrogen (H2) cannot lower the gas temperature below 
- 200K (Bromm, Coppi, & Larson 2002; Abel, Bryan & 
Norman 2002). Massive, metal-free stars shine at their Ed- 
dington luminosity, Le cx; M, and have roughly constant 
effective (surface) temperatures of ~ lO^K and lifetimes of 
~ 3xl0^yr, independent of their mass (Bromm, Kudritzki, 
& Loeb 2001). This implies that the number of ionizing 
photons produced per baryon incorporated into these stars 
(~ 4.5 X 10"') was independent of their mass function, and 
larger by more than an order of magnitude than the yield 
of the observed Pop-II (low-mass, metal-rich) stars. 

Did the first generation of metal-free stars produce a 
sufficient total number of ionizing photons to reionize the 
universe? This question is difficult to answer theoreti- 
cally. The UV radiation from the massive stars can eas- 
ily dissociate the fragile intergalactic H2 molecules which 
are responsible for star formation in mini-halos with virial 
temperatures Tvir < 10^ *K (Haiman, Rees, & Loeb 1997; 
Haiman, Abel & Rees 2000; Ricotti et al. 2002). More- 
over they could evaporate the gas out of the shallow po- 
tential wells of the first CDM mini-halos (Barkana & Loeb 
1999; Nishi & Tashiro 2000). Finally the hydrodynamic 
energy released by a single supernova explosion is much 
larger than the gravitational binding energy of the gas in 
these halos (Barkana & Loeb 2001). If feedback quenches 
star formation in mini-halos, then Pop-Ill star formation 
will be dominated by more massive halos where H2 can 
form in dense gas shielded from the UV background (Oh 
& Haiman 2002; Mackey, Bromm & Hernquist 2002). 

Analysis of the first year data from the WMAP satel- 
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lite suggests a large optical depth to electron scattering 
Tos — 0.17 ± 0.04, implying significant reionization at red- 
shifts as high as ~ 17 ± 5 (Kogut et al. 2003; Spergel et 
al. 2003). While this result is not surprising in view of 
earlier detailed calculations of the reionization history by 
metal- free stars (Wyithe & Loeb 2003; Paper I), it offers 
an intriguing empirical path for answering the aforemen- 
tioned question. In this Letter^ we argue that metal-free 
stars are required by the WMAP measurement of t^s- 

We assume the most recent cosmological parameters ob- 
tained through fits to WMAP data (Spergel et al. 2003). 
These include density parameters of f2m — 0.27 in matter, 
rjfc — 0.044 in baryons, JIa = 0.73 in a cosmological con- 
stant, a Hubble constant of Hq = 71 kms~^Mpc^^, an 
rms amplitude of erg = 0.84 for mass density fluctuations 
in a sphere of radius 8/i^^Mpc, and a primordial power- 
spectrum with a power-law index n — 1. Any tilt in the 
power spectrum at the level that is marginally indicated by 
the combination of the WMAP data with other data sets 
(Spergel et al. 2003) would reduce the small-scale power, 
delay reionization, and strengthen our conclusion. 

2. OPTICAL DEPTH TO ELECTRON SCATTERING 

Different reionization histories result in different large- 
scale polarization anisotropies of the microwave back- 
ground (Zaldarriaga 1997; Kaplinghat et al. 2002). 
A simple measure of the reionization history is the 
integrated optical depth to electron scattering Tcs = 

f^^ dz\cdt/dz\aTnc{z), where ctt is the Thomson cross- 
section and {dt/dz)'^ = —Hq{1 + z)\/Slm(l -I- z)^ -|- VL\. 
The electron number density, rie, depends on the mass fill- 
ing factors Q™ of the different ionized species in the gas, 

ne = (1 -f zf [Q-+n° + Q^,+ n°e + 2Q^^,^nl,] , (1) 

where vP denotes the total present-day number densities 
of hydrogen (H), and helium (He). 

3. REIONIZATION HISTORIES 

We now consider reionization histories under different 
model assumptions. The histories follow the reionization 
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of hydrogen as well as the single and double reionization 
of helium in an inhomogeneous IGM (Miralda-Escude et 
al. 2000). We include helium in our calculations since the 
harder spectrum of massive Pop-III stars can reionize both 
hydrogen and helium at high rcdshift (Oh et al. 2001; 
Venkatesan, Tumlinson & ShuU 2003). 

We calculate the emissivity of quasars described by the 
luminosity function of Wyithe & Loeb (2002a) as well as 
the emissivity of high redshift zero-metallicity stars and 
Pop-II stars (the latter is discussed in § 3.1). Our model 
has two free parameters: (i) the transition redshift, 2;tran, 
at which Pop-II stars start to dominate the cosmic radia- 
tion field of ionizing photons; and (ii) the product of the 
star formation efficiency in massive galaxies f^, (^stellar 
mass/total initial baryonic mass) and the escape fraction 
/esc of ionizing photons from galaxies. In Paper I we pre- 
sented contours of Tes as a function of 2tran and /*/esc- We 
have found that if the star formation efficiency in early low- 
mass galaxies resembled that in local dwarf galaxies with 
the same circular velocity, then values of /*/esc ^ 0.05 — 0.2 
and 2;tran ~ 15 lead to Tes within the range measured by 
WMAP. Below we summarize the basic ingredients of our 
model; more details can be found in Paper I. 

We treat the IGM as a two-phase medium. In regions 
of reionized IGM, star formation can only be initiated in 
halos that collapse above a virial temperature threshold 
of Tmin = 2.5 X lO^K, corresponding to the cosmological 
.Jeans mass of photo-ionizc;d IGM at ^ 10*K (e.g., Thoul & 
Weinberg 1996; Kitayama & Ikeuchi 2000). In the neutral 
IGM, star formation is limited to halos with virial tem- 
peratures above the minimum for efficient gas cooling. If 
only atomic hydrogen is present, then Tmin — lO^K. How- 
ever if molecular hydrogen is present then Tmin = 10^ '^K 
(Haiman, Thoul, & Loeb 1996; Haiman, Abel, & Rees 
2000) and star formation can proceed in much smaller ha- 
los. We begin by discussing the case with Tmin = lO^K. 

3.1. Ionizing Photons from Stars 

To compute the comoving densities of stellar ionizing 
photons (n-y) we integrate over model spectra for stellar 
populations to find the total number of ionizing photons 
(Nj) emitted per baryon incorporated into stars. We find 



larger fraction of their gas. Using the stellar mass TuUy- 
Fisher relation of Bell & De Jong (2001), we find the 
threshold circular velocity i)* = 176 km s~^ that corre- 
sponds to a stellar mass of 3 x W^^Mq at ^; = 0. The 
mass of a dark matter halo with circular velocity is 

^haio ~ 5 X 10^° [(1 + z) /10]~^/^ Mq. We define as the 
star formation efficiency in galaxies with circular veloci- 
ties larger than u*; below Vi, the star formation efficiency 
A (M/Mh*,ij'/'. Given T^i„ we find 
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pin is the mean comoving mass density of 
matter, and nst(> M) is the comoving number density of 
halos with masses > M (Sheth & Thormen 1999). Fi- 
nally, the value of dF^^/dz in the two phase IGM required 
for calculation of dn-^/dz may then be computed from 

dF, _ rff.(2.5 X lO^K) dF.jlO^K) 

Tz +(1-^H+) • (4) 

3.2. The Need for Metal-Free stars 

We begin with calculations of reionization histories in 
the absence of metal-free stars at all redshifts. We have 
calculated reionization histories for several values of /*/csc 
in terms of the volume filling factor Q of ionized regions 
in the low overdensity (< 20) IGM (Miralda-Escude et al. 
2000; Paper I), and computed the resulting values of Tes- 
Recently, Theuns et al. (2002) have argued that if HI 
reionization ended before z = 9, and the universe remained 
reionized after that, then adiabatic cooling would have 
lead to an IGM temperature at z = 3.5, much lower than 
inferred from the Lya forest. As a result, the most recent 
overlap of intergalactic HII regions must have occurred 
at 2; < 9, unless there was significant photo-heating at 
z 3.5 due to the double reionization of helium. In the 
absence of metal free stars we find that this constraint 
requires /*/esc ^ 0.1, resulting in a value of Tes < 0.09 
that is not consistent with the WMAP data. On the other 
hand, if we neglect the constraint that HII overlap oc- 
curred at z < 9, we find that in the absence of metal-free 
stars, reionization histories resulting in Tes ^ 0.13 imply 
HII overlap redshifts of ~ 15. However assuming feedback 
regulation of star formation in low mass galaxies, these 
histories require the unphysical value of /*/esc = 1.0 in 
galaxies with Vdr > u*. We therefore conclude that reion- 
ization by an early population of massive metal-free stars 
(with very large'^ iV^ ) must have contributed significantly 
to Tes at z ^ 15. 

The formation of metal free stars is moderated with cos- 
mic time by the gradual metal enrichment of the IGM. 
It has been argued that the dominant mode of star for- 
mation change to Pop-II once the metallicity of the as- 
sociated gas rises above IO^^Zq (Schneider et al. 2002; 
Bromm et al. 2001). If all Pop-III stars formed in the 
mass range of 14OM0 — 26OM0 where pair-unstable super- 
novae result in expulsion of metals into the surrounding 
gas (Heger & Woosley 2001), and if all metals were mixed 

^ We conservatively adopt the ionizing photon yield of IOOMq metal-free stars. The yield is larger by a factor of ~ 2 for Pop-III stars with 
masses > 300Mq (Bromm, Kudritzki, & Loeb 2001). The use of the higher yield would have lead to the same Tes for values of /*/esc at 
z > ztran that are smaller by a factor of 2 than we quote. 
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where rib is the number density of baryons, and i^*(z) is 
the fraction of baryons in the universe at redshift z that 
collapsed, cooled, and formed stars. 

A recent analysis of SDSS data (Kauffmann et al. 2002) 
has shown that in a large sample of local galaxies, the ra- 
tio between the total stellar {M^,) and dark matter halo 
(M) masses scales as M^/^ for < 3 x 10^°Mq, but 
is constant for larger stellar masses. Thus the star for- 
mation efficiency within dwarf galaxies is proportional 
to M"^/^. Since star formation is thought to be regu- 
lated by supernova feedback (Dekel & Silk 1986; Dekel 
& Woo 2002), and the supernova energy output is pro- 
portional to while the gravitational binding energy 
scales as Mti^;^ oc M^/-^, more massive galaxies retain a 
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uniformly throughout the entire IGM, then Pop-III star 
formation would have ceased as soon as ~ 0.3 ionizing 
photons were produced per cosmic hydrogen atom (Cia- 
rdi, Ferrara & White 2003). However, the highly inhomo- 
gcncous nature of the enrichment process (Furlanctto & 
Loeb 2003) and the formation of Pop-III stars beyond the 
above mass range, naturally lead to reionizion. 

Must the Pop-III stars possess a heavy initial mass func- 
tion (IMF)? TumMnson k ShuU (2000) computed the ion- 
izing photon output per baryon of zero-motallicity stars 
with a Salpeter IMF and found only a 50% increase over 
the values for metal enriched Pop-II stars with the same 
IMF. This moderate increase is to be compared with the 
increase by more than an order of magnitude for the 
massive Pop-III stars described by Bromm, Kudritzki, & 
Loeb (2001) (the number of ionizing photons produced per 
baryon is nearly independent of stellar mass). The argu- 
ments presented above therefore imply that Pop-III stars 
were very massive rather than of a Salpeter IMF. This 
provides the first observational support to the results from 
hydrodynamic simulations of early star formation (Bromm 
et al. 2002; Abel et al. 2002). 

In generating the star formation history we have extrap- 
olated results from the local galaxy sample, assuming that 
feedback governs the observed stellar mass as a function 
of halo potential depth. However high redshift galaxies 
may have cooling times that are shorter than the dynam- 
ical time of their cold gas. This may lead to a starburst 
that might render the time-delayed supernovae feedback 
less effective. We now briefly consider this case. Assum- 
ing the star formation efficiency to be independent of halo 
mass (case A in Paper I), we find that reasonable values 
of /*/esc ~ 0.1 can yield Tos ~ 0.17 in the absence of 
massive Pop-III stars. However, the corresponding reion- 
ization histories also have an early HII overlap. Therefore, 
unless the global product of escape fraction and star for- 
mation efficiency was more than an order of magnitude 
larger at ^; ~ 15 than at z ~ 9, our conclusion that an 
early generation of massive metal free stars is required to 
produce both t^s ~ 0.17 and HII overlap at ^; < 9 is in- 
dependent of whether feedback regulates star formation in 
low mass galaxies or not. 

3.3. Reionization Histories 

We now show that the joint constraints of t^s = 0.17 ± 
0.04 and HII overlap at z < 9 can be met simultaneously 
for values of 2;tran > 9 and /*/esc < 1 if massive metal- 
free stars were present at high redshifts. In Figure 1 we 
show two cases of the evolution of Q for hydrogen (solid 
lines) , singly ionized helium (dashed lines) and doubly ion- 



ized helium (dotted lines). We assume Ztran = 12 and 
/*/esc = 0.08 &i z < Ztran, yielding the most recent HII 
overlap at z ~ 9. A.i z > Ztran we assume /*/esc = 0.08 
(left) and /^/csc = 0.33 (right). The latter allows for the 
possibility of a higher star formation efficiency'' or escape 
fraction (Steidel et al. 2001) at high redshift. In the for- 
mer case there is a significant reionization but failed HII 
overlap at early times, followed by a recombination and a 
successful HII overlap at 2; ^ 9 due to Pop-II stars. The 

* Salvaterra & Ferrara (2003) have recently presented the observational case that 2;tran ^ 9 through fitting the near-infrared background. The 
lower value of ztran would allow a second HII overlap near z ~ 7 to be compatible with the large value of Tes and IGM temperature constraints. 

Note that the star formation efficiency in the low-mass galaxies of interest, r;* = 0.02/* (M/IO^M©)^/^ [(1 -|- z)/15], is generally < 



two examples yield Tes within one standard deviation and 
at the most likely value of WMAP, respectively. Note 
that our model leads to the double reionization of helium 
by quasars at z ~ 3.5 independent of the contribution 
from Pop-III stars. This is in excellent agreement with 
current observational data [e.g. Theuns et al. (2002) and 
references therein] . The reionization of Helll by quasars is 
characterized by a sharp rise in our model, with QHe++ in- 
creasing from ~ 50% at z ~ 4.3 to ~ 100% at z - 3.7. We 
do not expect pre- heating from the reionization of HcIII to 
impact on the conclusion that the most recent HII overlap 
took place at z < 9. 

Figure 1 also shows the evolution of Q for hydrogen 
when helium is ignored (light solid line). In these calcu- 
lations all photons blueward of 13.6eV reionize hydrogen. 
We find that the additional ionizing photons result in over- 
lap of hydrogen at an earlier epoch, but the net optical 
depth to electron scattering is decreased. 

So far we have assumed that atomic cooling dictates 
the mass of the smallest galaxies. However it has been 
suggested that significant production of X-rays could cat- 
alyze molecular hydrogen production allowing star forma- 
tion in mini-halos before reionization (Haiman, Rees, & 
Loeb 1996; Cen 2002). The light dashed lines in the left 
hand panels of Figure 1 show the evolution in Q for hydro- 
gen (in the absence of helium) assuming 

^min — 

lO^-^K 

in the neutral IGM. We find that the feedback-regulation 
built into the model sufficiently limits the star formation 
efficiency in mini-halos so that they do not contribute sig- 
nificantly to the reionization history. 

4. CONCLUSIONS 

We have shown that Pop-II stars alone cannot reion- 
ize the universe sufficiently early to produce the value of 
Tes ^ 17 ± 4% inferred from the WMAP data (Kogut et 
al. 2003; Spcrgcl et al. 2003). This conclusion holds if 
feedback regulates star formation in early low-mass galax- 
ies as observed in present-day dwarf galaxies. Without 
metal free stars, the condition that the most recent overlap 
of intergalactic HII regions must have occurred at z < 9 
(Theuns et al. 2002) generically yields t^s < 9%. The 
joint requirements that t^s ^ 17% and that HII regions 
overlapped at z < 9 imply that metal-free stars reionized 
a significant fraction of the IGM at z ~ 15. Moreover, we 
find that these stars must have been massive to provide 
the required order-of-magnitude increase in the production 
rate of ionizing photons relative to low mass Pop-II stars 
(Bromm, Kudritzki & Loeb 2001). Pop-III stars with a 
Salpeter IMF provide only a 50% increase in (Tumlin- 
son & ShuU 2000) and cannot account for the large t^s- 

For plausible reionization histories with metal-free stars 
we have found that singly and doubly ionized helium con- 
tributes as much as Ares = 0.02, and should be included 
in reionization calculations. We also found that if the star 
formation efficiency is rcgiilated in dwarf galaxies as ob- 
served in the local universe, then there is no significant 
contribution to the cosmic ionizing radiation field from ha- 
los with virial temperatures T^i^ < 10'''K, even if molecular 
hydrogen permits cooling of gas inside these mini-halos. 
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W.=0-08 (0-08) Zt„„=12.0 T„=0.13 W.=0-08 (0.33) Zt„„=12.0 T.,=0.17 




Fig. 1. — Sample rcionization histories assuming feedback regulated star formation, that produce large values of Tcb- Dark lines represent 
the evolution of the volume filling factors of ionized regions Qjj+' Qlle+ ^^'^ QHe++- ^1^° show calculations of Qjj+ in the absence of 
helium (light lines). The solid light lines show results for a pre-reionization minimum virial temperature for galaxies of Tmin = lO^K (atomic 
cooling) while the dashed light line corresponds to T^in = 10^ *K (molecular cooling). Each example is labeled by the values of /«/esc at 
z < ztran (and in parenthesis at z > ztran), ^tran and res- 



The evolution of the ionization state and metal enrieh- 
ment of the IGM during the reionization epoch can be 
best probed by infrared spectroscopy of gamma-ray burst 
(GRB) afterglows (Furlanetto & Loeb 2002). These tran- 
sient sources should be both abundant (Bromm & Loeb 
2002) and detectable (Lamb & Reichart 2000; Ciardi & 
Loeb 2000) out to high redshifts. The Swift satellite, 
planned for launch this year, is expected to detect ~ 100 
GRBs per year, out of which > 20% may originate at z > 5 
(Bromm & Loeb 2002). Photometric identification of the 
Lya trough can be used to select the high-redshift after- 
glows for which spectroscopic follow-up observations may 
be conducted on a 10-m class telescope (Loeb 2001). 

After the submission of this paper, a number of differ- 
ent groups have also commented on the implications of 
the large Tes measured by WMAP. Some studies find that 
Pop-II stars are able to reionize the universe sufficiently 
early to produce Tos > 0.13 (Ciardi, Ferrara & White 2003; 
Somerville & Livio 2003); however, this conclusion is chal- 
lenged by the existence of the observed Gunn-Peterson 



trough in quasars at z ^ 6.3 (as noted by Ciardi et 
al. 2003) and by the temperature evolution of the IGM 
(Theuns et al. 2002; Hui & Haiman 2003). Other authors 
agree that the combined requirements of a large Tes and a 
late HIT overlap implies that the efficiency of production of 
ionizing photons was higher at redshifts 2 > 10 (Haiman & 
Holder 2003; Sokasian, Abel, Hernquist & Springel 2003; 
Cen 2003). The conclusion that Pop-III stars needed to 
be massive to produce the observed t^s was also reached 
by Cen (2003) in agreement with our results. 

We thank Rennan Barkana, Volker Bromm, Steve 
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acknowledges generous support from the Institute for Ad- 
vanced Study at Princeton and the John Simon Guggen- 
heim Memorial Fellowship. This work was also supported 
in part by NSF grants AST-0071019, AST-0204514 and 
NASA grant ATP02-0004-0093. 
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